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ABSTRACT 

We find 9 nuclear cluster candidates in a sample of 14 edge-on, late-type galaxies observed with 
HST/ACS. These clusters have magnitudes (M/ ~ —11) and sizes {reff ~ 3 pc) similar to those found 
in previous studies of face-on, late-type spirals and dE galaxies. However, three of the nuclear clusters 
are significantly flattened and show evidence for multiple, coincident structural components. The 
elongations of these three clusters are aligned to within ~10° of the galaxies' major axes. Structurally, 
the flattened clusters are well fit by a combination of a spheroid and a disk or ring, with the disk 
preferred in two of three cases. The nuclear cluster disks/rings have F606W-F814W (^V-I) colors 
0.3-0.6 magnitudes bluer than the spheroid components, suggesting that the stars in these components 
have ages <1 Gyr. In NGC 4244, the nearest of the nuclear clusters, we further constrain the stellar 
populations via spectroscopy and multi-band photometry. This nuclear cluster is equally well fit 
by single-stellar populations with ages of either ^100 Myr or ~1 Gyr, and with masses of 2-5x10^ 
Mq. However, significantly better fits to the spectroscopy and photometry are obtained by combining 
two or more stellar populations. Exploiting emission lines that appear to originate ~1" from the 
NGC 4244 nucleus, we determine a lower limit on the dynamical mass of 2.bt\l X IO^Mq within 
19 pc, typical of values found for other nuclear clusters. We also present tentative evidence that 
another of the nuclear clusters (in NGC 4206) may also host a supermassive black hole. Based on our 
observational results we propose an in situ formation mechanism for nuclear clusters in which stars 
form episodically in compact nuclear disks, and then lose angular momentum or heat vertically to form 
an older spheroidal structure. We estimate the period between star formation episodes to be ~0.5 Gyr 
and discuss possible mechanisms for tranforming the disk-like components into spheroids. We also 
note the connection between our objects and massive globular clusters (e.g. lo Cen), ultra-compact 
dwarfs, and supermassive black holes. 

Subject headings: galaxiesmuclei - galaxies:star clusters - galaxies: spiral - galaxies:formation - galax- 
ies individual (IC 5052, NGC 4244, NGC 4206) 



1. INTRODUCTION 

Due to their deep gravitational potential wells and 
unique dynamics, the centers of galaxies frequently host 
exotic objects such as supermassive black holes, active 
galactic nuclei, and massive stellar clusters. Of these, 
only nuclear star clusters provide a visible record of the 
accretion of stars and gas into the nucleus. Therefore, 
studies of nuclear star clusters hold the promise to shed 
light on the formation of central massive objects of all 
types. 

Bright nuclear star clusters have been observed in 
^75% of the 77 face-on galaxies of type Scd and later 
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studied by iBoker et al.l l)2002|) . Similar clusters ap- 
pear to be frequent in earlier type spirals as well, 
although their study is complicated by the presence 
of bu lges and other nuclear phenomena ( CaroUo et all 
l2002j) . Dwarf elliptical galaxies, w hich may originate 
from low-mass spiral ga laxies (e.g. iMastropietro et all 
120051 iBeaslev et alJ2006D. also f requently po ssess nuclear 
star clusters l|Grant et alJl2005t lC6te et alJl2006 D. 

In late-type spiral galaxies with little or no bulge, nu- 
clear star clusters are very prominent and can be stud- 
ied in detail. From their s a mple of face-on, late-type 
spiral galaxies, IBoker et al.l l)2004() found that nuclear 
clusters have dimensions similar to Milky Way globular 
clusters, with typical half-light radii of 3.5 pc. How- 
ever, these clusters are significantly brighter than globu- 
lar cluste rs, with / band ma gnitudes ranging between —8 
and —16 ijBoker et al.ll200^ . Their high luminosities are 
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at least in part due to their large masses; iWalcher et al.l 
(jlpOS) have used high-resolution spectra to d erive typical 
masse s '^SxIO^Mq for nine galaxies in the iBoker et al.l 
l)2002() sample, simila r in mass to the most ma ssive Milky 
Way cluster, uj Cen l|Mevlan fc Mavo r '19861 However, 
the brightness of nuclear clusters also results in part from 
the presence of young stellar populations. Photometry 
and spectra of a number of clusters suggest that nu- 
clear clusters in both spiral and dwarf ellipticals have 
popula tions m uch youn ger than globular clusters (e.g. 

iHo et al. 1995t iLotz et al 2004^) . In late- type galaxies, 

most clusters appear to have stars with ages <100 Myr 
ijWalcher et al.l200fi^ . Furthermore, recent spectral stud- 
ies have shown that nuclear clusters are made up of com- 
posite stellar populations, with most having substan- 
tial old (>1 Gyr) stellar comp onent ( Lo ng et al. 2002; 
ISarzi et al.ll2005HWalcher et al.ll2 006: Ro ssa et al.ll200fl) . 

Not only are nuclear clusters present in a large frac- 
tion of present day spirals and dwarf ellipticals, they also 
have been invoked as the progenitors of massive globu- 
lar clusters and ultra-compact dwarfs. In a hierarchical 
merging scenario, many galaxies will be tidally stripped 
during merge rs leaving be hind only their dense nuclei 
(e.g. [Bekki & Chibal 1200^ . Nuclear cluster properties 
match the half-light radii and exceptional luminosity of 
massive globulars and UCDs ( M ackev & van den Bcrg^ 
11005; Philhpps et al. 2001; Jon es et all 12006). Further- 
more, the mul tiple stellar populations seen in uj Cen and 
Gl (e.g. Bedi n et al.ll20f)l lMevla~ al. 200l|, and the 
evidence for possible self-enrichment among the bright- 
est globular clusters in ellipticals (Harris et al. 2006) 
strengthen the association between these bright compact 
objects. 

Finally nuclear clusters appear to be intimately con- 
nected to the formation of supermassive black holes 
(SMBHs) at the centers of galaxies. For both 
spiral and elliptical galaxies, a clear correlation is 
seen between iiuclear cluster mass and bulge mass 
(iCote et al.l2006t|Wehner fc Harrisl2006HFerrarese et all 
l2006HR,ossa et al.ll2n06D . This correlation appears to fall 
along the well-known M, — a relatio n between black hole 
mass and bulge or galaxy mass fe.g. lFerrarese fc Merritti 
I200nt lFerraresdl2002.'l . Both SMBHs and nuclear clusters 
appear to contain the same fraction of the total galaxy 
mass, supporting the idea that they are two relate d 
types of central massive object (Wehncr fc Harrisll200^ . 
Whether a central massive object is a nuclear cluster or 
a SMBH appears to depend on mass; early-type galax- 
ies more massive than -10^° Mq appear to have SMBHs, 
while lower mass systems commonly host nuclear clusters 
(p^errarcse ct al. 2006). These findings strongly suggest 
a similar formation mechanism for nuclear clusters and 
SMBHs. Increasing evidence also points to massive star 
clusters as sites of intermediate-mass black hole (IMBH) 
formation ic.s . iMiller fc Co lbert"2004': 'Geb hardt et all 
H)05; Patruno et alJl200d). T hese IMBHs may also fol- 
low the Af, — a relation l)Gebhardt et al.. 2002). 

In this context, a detailed understanding of how nu- 
clear clusters formed may reveal the general mechanism 
driving the formation of central massive objects of all 
types. Two basi c scenarios have been suggested (e.g. 
IBoker et al.ll200^ : (1) nuclear clusters form from multi- 
ple globular clusters accreted via dynam ical friction (e.g. 
iTremaine et al.lll975t FLotz et al.ll200H) . and (2) nuclear 



clusters for m in situ from gas channeled into the center 
of galaxies l|Milosavlievicl2003) . In this paper we present 
observational evidence that favors the in situ formation 
scenario. In particular, we find nuclear clusters that are 
elongated along the major axis of the galaxies and have 
two coincident components. These clusters appear to 
consist of a younger flattened disk or ring and an older 
spheroid. Based on the evidence we present, we suggest 
a detailed formation mechanism in which nuclear clus- 
ters form in repeated star formation events from very 
compact nuclear gas disks. 

In 32 we describe the imaging and spectroscopic ob- 
servations used in the paper. Sj^l focuses on an analysis 
of the morphology and luminosities of the clusters as re- 
vealed by HST observations. The stellar populations of 
the multiple-component clusters are examined in 21 In 
[iniwe derive a dynamical mass for the NGC 4244 nuclear 
cluster, and then in i|Slpresent tentative evidence that the 
NGC 4206 nucleus may harbor both a star cluster and 
a supermassive black hole. We present a detailed forma- 
tion mechanism and discuss it in the broader context in 
IZI and then summarize our findings in fJS] 

2. OBSERVATIONS 

2.1. Identification of Nuclear Clusters 

Most of the data presented in this paper were taken 
with the Hubble Space Telescope's (HST) Advanced 
Camera for Surveys (ACS) as part of a Cycle 13 snap- 
shot program. Sixteen galaxies were observed in both 
the F606W and F814W filters, with exposure times of 
^700 seconds. Details of these observations and their 
reduction are presented in lSeth e t al. (2005) (Paper I). 

We identified nuclear clusters (NCs) by careful exam- 
ination of the central regions of our galaxies. The exact 
locations of the galaxy centers were determined using 
fits to 2MASS data presented in Paper I. Nuclear clus- 
ters were distinguished based on their proximity to the 
photocenter of the galaxy, their brightness/prominence, 
and their extent (i.e. if they appeared resolved). At the 
distance of most of our galaxies, NCs such as those ob- 
served by Bokcr ct al. ( 2002) should be at least partially 
resolved. 

Table lists the galaxies in our sample and whether 
each possesses a NC (shown in the fifth column). In 
two of the sixteen galaxies in the sample, the central re- 
gion was not covered by the ACS images (denoted with 
a "N/A" in Table Of the remaining 14 galaxies, 
seven have prominent and/or extended sources located 
within ~2 arcseconds of the photometric center of the 
galaxy. Another two have possible detections of NCs, 
where the candidate clusters were not very prominent 
(denoted with a "Maybe" ) . The remaining five galaxies 
have no obvious candidate NCs. Of these five definite 
non-detections, three (IRAS 06070-6147, IRAS 07568- 
4942, & NGC 4631) are in galaxies where strong dust 
lanes could have easily obscured any extant cluster. The 
other two of the non-detections (NGC 55 and IC 2233) 
are in low mass galaxies without strong dust lanes, sug- 
gesting that these galaxies truly lack NCs. 

In summary, we have detected candidate NCs in 50% to 
65% of our sample. Eliminating the galaxies with strong 
dust lanes, the fraction of galaxies with detected NCs 
may be as high as 82%. These detection fractions are 
consistent with other studies, both for late-type galax- 
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Fig. 1. — Images in the F814W filter of the nine candidate nuclear clusters. Each image is 2" on a side, and the bar in the bottom right 
corner indicates a length of 10 pc. Images have been rotated so that the x-axis lies along the position angle of the major axis of the galaxy 
disk. 



ies, where lB5ker et all l|2002|) finds ~75% of systems have 
clusters, and for early- type gal axies observed as part of 
the ACS Virgo Cluster Survey l)C6te et al.ll2006(l . How- 
ever, our sample is less well-suited than these studies for 
constraining the frequency of nuclear clusters because the 
edge-on perspective can result in more chance superpo- 
sitions and greater obscuration from dust. 

2.2. Apache Point Observatory 3.5-m Spectrum 

To better constrain the stellar populations of the 
NGC 4244 nuclear cluster (the nearest in our sample), 
we obtained long-slit spectra on May 19, 2005 with 
the Double Imaging Spectrograph (DIS) on the Apache 
Point Observatory (APO) 3.5m telescope. The DIS spec- 
trograph obtains both blue and red spectra simulta- 
neously. We used the "medium-resolution" blue grat- 
ing (1.22 A/pixel) and the "high-resolution" red grating 
(0.84A/pixel). The resulting images had good quality 



spectral coverage from 3700A to 5500A on the blue side 
and 5600A to 7200A on the red side. A 1.5" slit was 
used. Two exposures of 15 and 30 minutes were ob- 
tained centered on the nuclear cluster and oriented along 
the major axis of NGC 4244 with a position angle of - 
43° . We note that the elongation of the nuclear cluster 
is aligned to within ^10° of the galaxy's major axis (see 
Fig. Seeing at the time of observation was ~1.0", 
and thus the cluster was essentially unresolved. We re- 
duced these spectra using a PyRAF pipeline designed for 
DIS data reduction (written by Kevin Covey) which uti- 
hzes standard IRAF ONEDSPEC and twodspec routines. 
The NGC 4244 cluster observations were wavelength cal- 
ibrated using lamp spectra, and a further correction was 
applied to remove small shifts visible in the sky lines. 
The velocities were calibrated to the local standard of 
rest using rvcorrect and dopcor. Observations of 
Feige 66 at nearly identical airmass were obtained for 
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TABLE 1 

Galaxy Sample and Nuclear Cluster Properties 



Galaxy 


Type 


M^j Vmax 

[km/soc] 


Possess 
NC 


Multiple 
Components 




™F6061V - mp8i4w m-M 


IC 2233 


SBcd 


-16.60 


79 


No 


N/A 






29.96 


IC 5052 


SBcd 


-17.23 


79 


Yes 


Yes 


-10.77 


1.08 


28.90 


IRAS 06070-6147 


Sc 


-20.50 


133 


No;Dust 


N/A 






31. 22^^ 


IRAS 07568-4942 


SBc 


-17.70 


142 


No;Dust 


N/A 






30.06<: 


IRAS 09312-3248 


Sc 


-17.90 


98 


Maybe 


No 


-9.24 


1.36 


30.61^^ 


NGC 55 


SBm 


-17.05 


67 


No 


N/A 






26.63=^ 


NGC 891 


Sb 


-18.78 


214 


N/A 


N/A 






29.61*= 


NGC 3501 


Sc 


-17.37 


136 


Yes 


No 


-10.58 


1.62 


31. 04^ 


NGC 4144 


SBc 


-17.20 


67 


Maybe 


No 


-8.33 


0.85 


29.36=^ 


NGC 4183 


Sc 


-18.41 


103 


Yes 


No 


-11.15 


0.93 


31.35* 


NGC 4206 


Sbc 


-18.46 


124 


Yes 


Yes 


-14.17 


0.93 


31. 29s 


NGC 4244 


Sc 


-17.08 


93 


Yes 


Yes 


<-12.65' 


0.77 


28.20=^ 


NGC 4517 


Sc 


-18.18 


137 


Yes 


Maybe 


-11.34 


1.77 


29. 30*= 


NGC 4565 


Sb 


-20.47 


227 


N/A 


N/A 






31. 05^ 


NGC 4631 


SBcd 


-19.75 


131 


No;Dust 


N/A 






29.42=^ 


NGC 5023 


Sc 


-16.29 


77 


Yes 


Maybe 


-9.79 


1.26 


29.10=^ 



Note. — Types and velocities from HYPE RLEDA/LEDA iPaturel et al.lll93Rl I20n,'jl^. ( a.) | ^eth et al.l 120051) (b) 
IBottinelli et al. I 1985), fc) "Bottin olli et all 119881. ( d) ass umins Ho=70 km/sec, fe) ITonrv et al.l i200lt . ffl TOllvfc Pierce! 
i200Cll . (gTlGavazzi ct al. 1 1999). (h) IJensen et al.l 1200311 (i) The derived magnitudes on the NGC 4244 cluster are lower 
limits due to the central pixels being saturated in both bands. 



flux calibration shortly before the observations; however 
the night was not photometric and thus the absolute flux 
scale is not fully constrained. The two observations of 
the NGC 4244 cluster differ in flux scale by under 10%, 
and our absolute flux is probably uncertain by a simi- 
lar amount. Combined, the two spectra give a S/N of 
^^35 per pixel. The cluster spectrum was extracted from 
within a 2" diameter aperture, and regions used to de- 
fine the background subtraction were within the confines 
of the galaxy, 20"-30" (-l/4th the scale length of the 
galaxy) along the major axis in both directions. Any 
contamination from the underlying galaxy component 
should be very small, as the source was well above the 
background at all wavelengths. We analyze this data in 
§4 and §5. 

3. NUCLEAR CLUSTER MORPHOLOGIES AND 
LUMINOSITIES 

The most basic result of this paper is that nuclear 
clusters are not all simple, single-component objects, 
as is evident from inspection of Figure ^ which shows 
the F814W images of all nine nuclear cluster candi- 
dates. Of these candidates three (IC 5052, NGC 4206, 
&; NGC 4244) look like miniature SO galaxies, possess- 
ing both an elongated disk or ring component and a 
spheroidal component. These components are both com- 
pact with physical dimensions of tens of parsecs or less. 
We will show in f|31l.3 that these components appear to 
physically overlap. We therefore refer to these clusters 
as multi-component clusters throughout the rest of the 
paper. Two other galaxies, NGC 4517 & NGC 5023 also 
show tentative evidence for similar structures. In this 
section we focus on modeling the morphology of the clus- 
ters, in particular those that have multiple components. 
We will show: (1) the nuclear cluster elongations are 
closely aligned with the disk of the host galaxy, (2) mod- 
els with both a spheroidal and disk or ring components 
fit the multi-component clusters significantly better than 
single-component models, and (3) the nuclear clusters in 



our sample have magnitudes and sizes similar to those 
observed in other galaxies. 

3.1. Model fitting 

The dimensions and scales of the clusters were deter- 
mined via 2-D fitting to analytical functions (e.g. King 
profiles). Because the objects are very compact, it was 
necessary to include the effects of the PSF. We used two 
similar methods to fit the nuclear clusters. First, we used 
the i shape program ( L arsen . 1 999i) , which was designed 
to fit semi-resolved clusters. Unfortunately, ishape only 
works with elliptically symmetric objects, which some of 
our objects clearly are not. Therefore, we wrote our own 
code in IDL, using the existing MPFIT2DFUN package to 
perform Levenberg-Marquardt least squares fits. 

In both ishape and our program, analytical functions 
are convolved with the PSF before fitting to the data. 
We performed this convolution using a position-variable 
PSF derived from our data (see Paper I for details), and 
over-sampled by a factor of 10 to minimize the effects 
of aliasing and binning. A separate PSF was generated 
for each cluster based on its position on the chip. We 
conducted tests to insure that our program give similar 
results to the ishape program. 

We used ishape (1) to determine what types of func- 
tions best fit the data, and (2) to assess if each clus- 
ter is indeed resolved. We fit Gaussian, King, Moffat, 
and Hubble profiles to each cluster and found that King 
and Moffat profiles provided the best overall fits in terms 
of the reduced x^, as also found by Boker ct al. (2004). 
Both profiles provided equally good fits to our clusters, 
and we chose to continue our fitting with elliptical King 
profiles of the form: 

E(z) = So|^^^-^^^=| (1) 
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where is the projected surface brightness, x and 

z are the coordinates along the major and minor axes, 
and q is the axial ratio (minor/major). We have fixed 
the concentration (c = rudailTcore) to be 15 for our 
fits, based on the fit concentrations of our brightest clus- 
ters and simil ar to that found for Galactic analog uo Cen 
ijTrager et alJ ll995). We report only half-hght/effective 
radii as this quantity is well det ermined even for rel- 
ative ly low signal-to-noise data l)Carlson fc HoltzmanI 
120011) . and docs not depend strongly on the chosen model 
(see Fig. 2, Bokcr ct al. 2004). 

Given sufficient signal-to-noise, semi-extended sources 
can be d etected with i ntrinsic widths as small as 10% of 
the PSF l|Larsenlll999t) . The FWHM of our PSF is -0.9 
pixels, meaning objects with widths smaller than a pixel 
should be easily resolvable. The signal-to-noisc of our fits 
ranged from '^50-1500. For four of the galaxies (IC 5052, 
NGC 4206, NGC 4244, & NGC 5023), the FWHM of the 
fitted ishape King profiles is greater than one pixel, and 
the objects are clearly partially resolved. The values 
of these fits is 2-20 times better than the of a point- 
source fit (see last column, Table 0). The remaining fits 
have FWHM >0.2 pixels, and have values 1-3 x the 
best-fitting point-source . Of these, NGC 4144 and 
NGC 4517 are stiU clearly resolved, while IRAS 09312- 
3248, NGC 3501, & NGC 4183 (all of which are at large 
distances) are only marginally resolved. 

3.1.1. Elliptical King Model Fits 

Table|21shows the results of our program's fits to single- 
component elliptical King profiles with c = 15. These fits 
show: 

1. All the well resolved clusters are flattened, as 
shown in Figure[21 The median axis ratio [q = h/ a) 
is 0.81, with q - 0.4 for NGC 4206 and NGC 4244, 
the two systems with the most prominent disk (see 
Sll.2). 

2. The position angle of the three multi-component 
nuclear clusters (IC 5052, NGC 4206, & NGC 4244) 
are all aligned within 10 degrees of the galaxy disk 
position angle. This can be seen plainly in Fig- 
ures [Hand □ The NGC 4517 and NGC 5023 clus- 
ters also are somewhat elongated along the plane 
of the galaxy. 

3. The effective radius (also known as the half-light 
radius) ranges between 1 and 20 pc, with most of 
the clusters having effective radii between 1-4 pc. 
Figure |3I demonstrates that the measured effective 
radii are similar to those observed for other ob- 
jects hi£ludiiigjate;_t^ nuclear clusters (filled cir- 
cles [BSte^FaO |200^_aiid_galactic globular clus- 
ters (open triangles lHarrislll996 D . 

4. After subtraction of the King profile fit, the resid- 
uals of IC 5052, NGC 4206, and NGC 4244 clearly 
show evidence for a flat disk- or ring-like compo- 
nent (see the middle column of Fig. 0)). The re- 
duced of these three galaxies is also very high. 




0.2 0.4 0.6 0.8 

q (minor/major) 



1.0 



Fig. 2. — Elliptical King profile fit results: the axial ratio q is 
plotted against the position angle of the major axis of the nuclear 
cluster relative to the galaxy disk (APA). The connected data 
points connect the fits for each galaxy in the two filters. Open 
circles indicate the three galaxies (IRAS09312-3248, NGC 3501, 
NGC 4183) for which good fits were not achieved, due to the NC 
compactness and complexity of the surrounding emission. This 
figure shows that most of the clusters we have detected were 
significantly flattened, and that the flattening aligns closely with 
the plane of the galaxy. 



suggesting that the single component King pro- 
file is not a good fit to the data. We note that 
NGC 4517, which has a disk-like appearance in 
its outer isophotes, also has a high x^ value. In 
NGC 4244 the central eight pixels were saturated 
and therefore excluded from the fit. 



3.1.2. Elliptical King + Disk Models 

To account for the fiattened residuals remaining af- 
ter fits to the elliptical king profiles in the IC 5052, 
NGC 4206, and NGC 4244 nuclear clusters, we have fit 
an additional morphological component to these clusters. 
In this section we examine fits to disk models, in the next 
section we consider fits to ring models and compare the 
disk and ring fits. 

As a disk model we adopted the surface brightness pro- 
file o f a transparent, self-gravitating, edge-on exponential 
disk l)van der Kruit fc Searle 1981) : 



E(x, z) = T,o—Ki 



sech 



(2) 



where x and z are the major and minor axis coordinates, 
Eq is the central surface brightness, Ki is the modified 
Bessel function, is the disk's exponential scale length 
and zq is the scale height of the disk. We note that we 
tried other vertical profile fits (e.g. simple exponential) 
which were too concentrated near the midplane to pro- 
vide good fits to the data. 

Figure 01 shows that the inclusion of a disk compo- 
nent clearly improves the fit relative to the single com- 
ponent elliptical King profile. Quantitatively, the x^ of 
the fits shown in Table 13 is significantly reduced relative 
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TABLE 2 
Elliptical King Profile Fits 



Galaxy 


Pix. Scale 


Filter 


' 1 


n 
H 


APA^ 


Fteduced 


A.PS/ A 










(0/a) 




X 




IC 5052 


1.46 


F606W 


3.64 


0.76 


-3.8 


23.8 


5.78 


IC 5052 


1.46 


F814W 


2.93 


0.80 


-6.6 


18.1 


7.00 


IRAS 09312-3248 


3.21 


F606W 


2.46 


0.68 


43.1 


2.8 


1.50 


IRAS 09312-3248 


3.21 


F814W 


2.24 


0.86 


39.8 


3.5 


1.42 


NGC 3501 


3.91 


F606W 


6.82 


1.00 


89.9 


8.2 


1.41 


NGC 3501 


3.91 


F814W 


3.58 


0.61 


19.2 


10.1 


1.38 


NGC 4144 


1.81 


F606W 


1.46 


0.85 


83.4 


3.2 


2.69 


NGC 4144 


1.81 


F814W 


1.99 


0.81 


74.8 


4.6 


2.30 


NGC 4183 


4.51 


F606W 


2.76 


1.00 


-40.8 


9.6 


1.20 


NGC 4183 


4.51 


F814W 


3.02 


0.78 


68.8 


6.4 


1.12 


NGC 4206 


4.39 


F606W 


16.30 


0.37 


5.6 


54.7 


2.39 


NGC 4206 


4.39 


F814W 


15.81 


0.40 


5.3 


42.4 


2.75 


NGC 4244 


1.06 


F606W 


2.92 


0.46 


-10.0 


94.1 


20.73 


NGC 4244 


1.06 


F814W 


3.35 


0.49 


-9.9 


82.3 


19.63 


NGC 4517 


1.76 


F606W 


1.71 


0.85 


-22.9 


46.0 


1.20 


NGC 4517 


1.76 


F814W 


1.33 


0.85 


-16.1 


64.2 


1.23 


NGC 5023 


1.60 


F606W 


12.54 


0.81 


16.2 


13.4 


4.27 


NGC 5023 


1.60 


F814W 


10.22 


0.83 


18.6 


11.1 


6.89 



Note. — (a) APA is the position angle relative to the galaxy disk. 
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Fig. 3. — The fitted effective radius vs. the absolute aperture magnitude in the F814W filter. Large filled and open circles as in Fig. 1^ 
Small symbols de note the following: Filled circles - late-type nuclear clusters from[B6kcr ct al] 120041) . Open circles - dE,N nuclei from 
IDe Propris et aTl (20050 (transformed to Mj using Padova models which suggests stellar populations older than 1 Gyr have F814W-F850LP 
colors of 0.25it0.1 mags). Filled tr iangles - ma ssive globular clusters in NGC 5128 observed by Martini & Ho (2004). Open triangles - 
Galactic globular clusters from the IHarnl 119961) catalog (note that ui Cen is the brightest of these). Filled Stars - Forna x UCDs from 
l^inkwatcr ct al. ( 2003), assuming V-I of l.litO.l as shown in Micskc et al. 1 2002). Open Stars - faint fuzzy M31 clusters from lHuxor et all 
I2b05i) . Open Squares - young massive clusters (YMCs) from iBastian et al.i t2009) . Note that Padova models suggest Mpgi^w ~ Mj to 
within a few hundreths of a magnitude. 
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Fig. 4. — Residuals of model fits to F606W observations of 
nuclear clusters in IC 5052, NGC 4206, and NGC 4244. The 
images show the area over which fits were done, and are in the 
coordinate system of the ACS images. They are thus rotated 
from the images in Figs. [Tl and ISl Bars in the upper right show a 
length of 10 pc in each galaxy. The left most column shows the 
original image, the middle and right columns show the images 
after subtraction of the best-fitting elliptical king and elliptical 
King + disk profiles. All are shown with the same scaling. In 
each case, the elliptical King -|- disk profile clearly provides a bet- 
ter fit to the cluster than a single-component elliptical King profile. 



to those for the single-component fit (Table 12) • The im- 
provement is most dramatic in NGC 4244, and is least 
significant in IC 5052. However, the reduced values 
remain well above unity suggesting the model doesn't fit 
the data perfectly. For IC 5052 and NGC 4244, the disk 
scale lengths are comparable to the effective radii of the 
spheroidal components {hx ^ 3 pc), while in NGC 4206 
the disk component is very extended (hx = 31 pc). 

In NGC 4206 and IC 5052, the residuals show a similar 
pattern which suggests the presence of a compact central 
source and/or a ring rather than a disk. We discuss the 
issue of whether the components are rings or disks in 
the section below. Both these galaxies also show a color- 
gradient perpendicular to the major axis, suggesting the 
presence of a thin dust disk. 

3.1.3. Rings vs. Disks 

In this section we present evidence that the flattened 
components in NGC 4206 and NGC 4244 clusters are 
better fit by a disk model than a ring model, while 
in IC 5052, a ring model provides a better fit. In all 
three cases, it appears that the disk or ring extends to 
small radii and overlaps with the light distribution of the 
spheroidal component. 

In Tabled we present the best- fitting ring + elliptical 
king models for each image. For the ring component, 
we assumed a constant midplane density with an outer 
(fout) and inner (ri„) radius: 

if X < : (3) 



E(x,z) = 2posech' (j-^ (^^JrLt-^'^ 
else : ^{x, z) = 

This model is brightest for Vm < x < rout, and con- 
tributes less flux near the center, making it qualitatively 
different from the disk model (Eq. 2). However, given 
that the details of both models are somewhat arbitrary, 
the results of the fits should be viewed with some caution. 

For the NGC 4206 and NGC 4244 NCs the ring fits 
(Table^J have reduced values that are 20-75% higher 
than the disk values shown in Table 13 Furthermore, the 
best-fitting ring models have inner radii extending all the 
way to the center of cluster. Thus, the evidence suggests 
that in the NGC 4206 and NGC 4244 NCs, the second, 
flattened component is disk-like rather than being ring- 
like. 

In contrast, the NC in IC 5052 shows two symmetric 
spots in its residuals (Fig.0J and color map (Fig.EJ. Not 
surprisingle, the ring model provided a better fit (20% 
decrease in reduced x^)i with an inner and outer radius 
of the ring of 6.8 pc and 11.8 pc, respectively. The whole 
ring falls well within the tidal radius of the spheroidal 
component (~24 pc), and thus their stellar distributions 
appear to overlap. 

Nuclear ring s have previo usly b een d etected in molec- 
ular gas (e.g. [B okcr ct al. 1997, Il998j) and UV stellar 
light (e.g. lMaoz et al. 1996') . However, these studies find 
rings an order of magnitude larger in size than the one we 
find in IC 5052. The only previously observed structure 
around a nuclear cluster on a similar scale to the IC 5052 
ring and NGC 4244 disk is the compact molecular disk 
detected coincident wi th the nuclear cluster in IC 342 
ijSchinnerer et alJl2003j) which has a diameter of 30-55 pc 
depending on the assumed distance to the galaxy. We 
note that this disk is barely resolved by their CO obser- 
vations (with a resolution of 1.2"), and may therefore be 
either a disk or a ring. 

The evidence in all three clusters that the disk or ring 
components do overlap with the spheroidal components 
justifies our labeling these three nuclear clusters as multi- 
component clusters. For simplicity, and because disk 
models are prefered in two of the three multi-component 
clusters, we will refer to the flattened components in the 
next few sections as being disks. We will return to the 
formation and possible role of rings in our discussion 
(331.3). 

3.2. Luminosities 

The F814W absolute magnitude {Mpsiiw) and the 
F606W-F814W color are given in Table U] Photometry 
was done using t he Vega-based z eropo ints and reddening 
coefficients from iSirianni et all l)2005(l . The foreground 
reddenings we r e corre cted using values of E{B — V) from 
iSchleeel et all fm^), which were at most 0.05 mags. 
These measurements are complicated by the high and 
variable background level and the complex shapes of the 
nuclear clusters. A scatter of ~0.1 magnitudes was found 
by varying the aperture size within reasonable limits. 
The formal errors on the calculated magnitudes arc quite 
small (typically <0.01 mags). Magnitudes from the fit- 
ted models agree with the aperture magnitudes in cases 
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TABLE 3 

Elliptical King + Disk Profile Fits 



Galaxy 


Pix. Scale 


Filter 




Disk 




Ell. 


King 




APA 


Reduced 




[pc/pixcl] 




Mo 


hx 




fj-o 




<? 


[°] 










[mag/n"] 


[pc] 


[pc] 


[mag/n"] 


m 








IC 5052 


1.46 


F606W 


17.27 


3.54 


1.47 


16.15 


3.97 


1.00 


-10.8 


20.3 


IC 5052 


1.46 


F814W 


16.28 


3.48 


1.19 


14.74 


3.55 


1.00 


-12.0 


15.7 


NGC 4206 


4.39 


F606W 


17.44 


31.00 


15.76 


12.87 


4.94 


0.45 


5.1 


26.8 


NGC 4206 


4.39 


F814W 


16.61 


31.36 


17.42 


12.00 


5.32 


0.44 


5.0 


16.7 


NGC 4244 


1.06 


F606W 


13.64 


2.51 


1.64 


16.29 


10.79 


0.78 


-10.1 


15.1 


NGC 4244 


1.06 


F814W 


13.18 


2.73 


1.42 


13.84 


5.70 


0.73 


-9.6 


16.3 



TABLE 4 

Elliptical King -t- Ring Profile Fits 



Galaxy 


Filter 




Ring 




Ell. 


King 




APA 


Reduced 














Mo 




<? 


[°] 








[pc] 


[pc] 


[pc] 


Ltot 


[mag/n"] 


[;if 








IC 5052 


F606W 


11.83 


6.83 


1.44 


0.18 


15.64 


3.24 


1.00 


-11.0 


16.9 


IC 5052 


F814W 


12.13 


6.91 


1.61 


0.14 


14.36 


3.05 


1.00 


-11.5 


13.2 


NGC 4206 


F606W 


84.82 


0.44 


18.56 


0.53 


13.57 


7.69 


0.48 


5.3 


32.5 


NGC 4206 


F814W 


83.78 


0.44 


20.69 


0.53 


12.59 


7.70 


0.47 


5.2 


21.5 


NGC 4244 


F606W 


10.42 


0.11 


0.92 


0.18 


13.19 


4.07 


0.58 


-10.1 


26.7 


NGC 4244 


F814W 


10.75 


0.14 


0.91 


0.16 


12.51 


4.25 


0.59 


-9.7 


23.5 



where the data is well- fit, but differ substantially when 
the fits are poor. 

Figure 131 plots the F814W aperture magnitudes mea- 
sured for all nine nuclear candidates against their fitted 
half-light radius from the elliptical King profiles. We 
also plot magnitudes and sizes of other compact stel- 
lar systems, assuming that Mj ~ MpsiAW since the 
F814W filter is very similar to the Johnson I filter, with 
single-stellar populations expected to be '^O.OS magni- 
tudes brighter in F814W than in I l|Girar di 2006). 

Examination of Figure |31 shows that the magnitudes 
and sizes of t he objects in our sample are similar to the 
NCs found bv lRoker et all ^M^k ifiHed circles) for 
their sample of face-on, late-type galaxies. |B6ker et aO 
l)2004tl measure sizes for only the brighter clusters in their 
sample, but find NCs with measured magnitudes as faint 
as Ml = —8.6. Because of their galaxies' face-on ori- 
entations, clusters with disk and spheroid components 
such as we have observed would have gone unnoticed in 
their study. The morphology of the M33 nuclear cluster 
is flattened (e = 0.17 ) along an axis parallel to the ma- 
jor ax is of the galaxy ijLauer et al.lll998t llVIatthews et al.l 
119991) ■ as would be expected for a multi-component clus- 
ter observed at the inclination of JVISS's disk. The stellar 
popu lations in the spectra of both ]VI33's NC (Long et alj 
I2OOI and som e of the objects in the .Boker et al.. (.200211 
sample ( Walch er et aLllToOffi . are also similar to what we 
see in igjfor NGC 4244's NC. It is therefore clear that 
the nuclear clusters observed here are similar to those 
found previously. Furthermore, the multi-component 
morphology may be a common feature of nuclear clus- 
ters. 

We also note that the two most luminous clusters 
are those with the most prominent disk components 
(NGC 4206 & NGC 4244). This correlation is not sur- 
prising, as the disks seem to be composed of recently 
formed stars (see 21 ^^'^ would therefore have low mass- 



to-light ratios. Over time, as the young stellar popula- 
tions age the disks will become less prominent. 

The observed color of the clusters (see Table ^1 can 
constrain the total dust column to the center of the host 
galaxies. We estimated rough lower and upper limits 
on the total reddening by comparing the colors of the 
NCs to that of an old and red, or young and blue stellar 
population. For most of the galaxies, the F606W-F814W 
NC colors of ~1 suggest that the reddening lies between 
E(B-V) of and 1. However, for NGC 3501 & NGC 4517 
the reddening must be greater than 0.7 mags. 

4. NUCLEAR CLUSTER STELLAR POPULATIONS 

In this section we demonstrate that the multiple struc- 
tural components seen in some NCs in 331 also have dis- 
tinct stellar populations. In particular, the disk com- 
ponents are made of younger stars than the spheroidal 
components. We first examine color maps of the clusters 
and argue that the disks most likely have populations 
younger than 1 Gyr. We then examine the bright nearby 
NC in NGC 4244 in detail, using multi-band photometry 
and spectroscopy (see Fig. 01 to show that the spectra is 
best fit by two or more epochs of star formation. 

Two recent papers have presented very thorough stud- 
ies of the stellar populations of nuc lear clusters in both 
early and late typ e spiral galaxies (| Walcher et al.ll200^ 
iRossa et alJ2006li . Both use fitting of spectra to convinc- 
ingly show that nuclear clusters contain multiple gen- 
erations of stars. For a sam ple of nine bright n uclear 
clusters in late-type galaxies. iWalcher et al.l f2006) has 
shown that all of the clusters observed had stellar popu- 
lations younger than 100 JVIyr. The typical mass of the 
youngest stars was a few x 10^ IM©. However, in most 
of the clusters, the bulk of the stellar mass had an age 
of >2.5 Gyr. Our findings on the stellar populations in 
our multi-component clusters, presented below, are fully 
consistent with these papers. 
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Fig. 5. — Maps of the F606W-F814W [VEGAMag] color for each cluster. Contours showing the F606W brightness are overlaid in black. 
The clusters with obvious disk components (IC 5052, NGC 4206, & NGC 4244) all show bluer disks and redder spheroids. The color bars 
on the right side indicate values for the F606W-F814W color. Colors have been corrected for foreground reddening. The central eight 
pixels of the NGC 4244 NC were saturated and thus don't give accurate colors. The circled region in NGC 4244 shows the suspected HII 
region used to obtain a dynamical mass estimate in JS] 



4.1. Color Maps 

FigureElshows F606W-F814W color maps for the three 
nuclear cluster candidates identified as having multiple 
components. In each cluster, an obvious and significant 
color difference exists between the disk and spheroidal 
components. The disks are 0.3-0.6 magnitudes bluer in 
F606W-F814W than the reddest parts of the spheroids. 

The observed color difference between the disk and 
spheroid components is most simply interpreted as a dif- 
ference in stellar population. Despite the unknown red- 
dening, we can place upper limits on the age of the disk 
compon ent. Using th e Padova single stellar population 
models ()Girardill2006() . for a 10 Gyr old stellar popula- 
tion with metallicities between [Fe/H] of -1.7 and -0.4, 
the F606W-F814W color ranges from 0.73 to 0.91, with 
redder colors for more metal-rich populations. The dif- 
ference in color between a 10 Gyr and 1 Gyr population 
is ^0.25 magnitudes regardless of metallicity, suggesting 
that disks cannot be composed solely of intermediate age 
(>1 Gyr) stars. At ages younger than 1 Gyr, however, 
the color quickly becomes bluer. Thus the > 0.3 mag 
color difference is evidence that the disks have stellar 
populations younger than 1 Gyr. Put another way, we 
find that the disk components all must have intrinsic col- 
ors bluer than F606W-F814W=0.65, which corresponds 
to single-stellar populations younger than 1 Gyr. Our 
analysis of the spectrum of NGC 4244's NC con- 
firms that stars with ages <100 Myr are present. Fur- 
thermore, each of the nine nuclear cluster spectra ana- 



lyzed bv lWalcher et al.l l)2006j) showed evidence for stellar 
populations with ages <100 Myr 

The suggested interpretation of the color difference as 
a difference in stellar population may be complicated by 
a number of factors. First, even if two separate compo- 
nents exist, both likely contribute to the observed color 
at most locations. This overlap suggests that the actual 
spread in age is larger than implied by the color differ- 
ences. Second, the colors could be affected by the pres- 
ence of dust structures on the same scale as the clusters 
or by emission line contribution. Based on the symme- 
try of the nuclear clusters (e.g. the disk is roughly the 
same color on both sides of the nucleus), and the spectral 
results from ti4.2l we believe that the color difference is 
unlikely to be due solely to dust or emission lines and 
instead really does reflect a difference in stellar popula- 
tion. 

4.2. Stellar Populations in NGC J^Hi 

In the previous section we showed that the color off- 
set between the disk and spheroid components of the 
nuclear cluster indicates that the disk has a systemati- 
cally younger stellar population. However because of the 
uncertain reddening, the ACS colors alone cannot accu- 
rately constrain the ages of the stellar populations. In 
this section we use our APO spectrum and photometry 
of the NC in NGC 4244 from the Sloan Digital Sky Sur- 
vey (SDSS) and 2MASS survey to constrain the cluster's 
age. Before describing the details of this analysis we first 
describe the photometry. 
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Fig. 6. — A portion of the APO/DIS spectrum of the NGC 4244 NC with best-fitting models overplotted. Only the portion of the 
spectrum with the most significant spectral features are shown. Top Panel - single-age best fits at 70 Myr and 0.8 Gyr - the ages at 
which there are minima in the joint-x'^ distribution (see Fig.|3- The most age-sensitive portion of the spectrum lies around the Ca H line 
at ~3940A, although numerous other features also cannot be fit well with a single-age population. Middle Panel - the best two-age fit, 
including both stars with ages of 0.1 Gyr (dark gray line) and 1.0 Gyr (light gray line). This fit is a significant improvement relative to 
either single-age fit. Bottom Panel - Residuals for all three fits. 



The nuclear cluster in NGC 4244 is prominent, and 
stands well above the emission from the galaxy disk. 
Because of this, it was detected as a point source in 
both the SDSS and 2MASS surveys. This data com- 
pliments the spectroscopic data by giving a large wave- 
length baseline to constrain th e mass, reddening, and 
ages of the NC dde Criis et a].l l2n03V In both surveys, 
the source is essentially unresolved, and photometry is 
determined from PSF-fitting. In the SDSS survey, the 
DR4 object ID is 58773909806344622, with PSF mag- 
nitudes of 17.81, 16.62, 16.18, 15.96, and 15.73 in the 
u', g', r', i', and z' filters. In 2MASS, the object ID is 
12172945-^3748264, with magnitudes of 14.36, 13.75, and 
13.55 in J, H, and K bands. The errors on the SDSS 
magnitudes are ~0.G2 mags, while the 2MASS observa- 
tions have errors of '^0.05 mags. The difference in res- 
olution and photometric methods between the two sur- 
veys means there may be some relative systematic error 
in the magnitudes. Therefore, any conclusions based on 
the photometry alone should be considered with caution. 
We note that because the central pixels of the NC were 



saturated in the ACS data (taken with gain=l e~), we 
were unable to derive accurate photometry from the HST 
images. 

4.2.1. Single-Age Fits 

In principle, the photometry and spectrum can be 
fit with arbitrary stellar populations, but in practice 
this leads to far too many degeneracies. To simplify 
the problem we started by assuming the NGC 4244 
cluster contains stars of a single age and then deter- 
mine a best-fitting extinction, mass and ^is a func- 
tion of age from both the photometric and spectroscopic 
data separately. For the photometric data, we used the 
single-stellar popula tion model s from the Pa dova group 
ijcirardi et aLll2004|) using the lMarigol l)200H) treatment 
of the AGB. These models have been produced for both 
t he SDSS and 2 MASS photometric systems and assume 
a lKroupal 1)200 Ij) initial mass function (IMF) between 0.1 
and 100 Mq. The extinction was cal culated in A„ assum- 
ing R„ = 3.1 using coefficients from Cirardi et al. (200^ 
and lCardelli et all l[l989 ) for the Sloan and 2MASS data 
respectively. For the spectroscopic fitting, single stellar 
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population spectral templates from iBruzual fc Chariot 
1 were used. These templates assume a Tchabriei 
1 2003<) IMF from 0.1 to 100 Mq. They have lA resolu- 
tion, and were redshifted and interpolated to the wave- 
lengths of our observed spectrum. The reddening was 
apphed to the spectrum using the iCardelli et al.l l|1989D 
prescription. We fit the feature rich blue side of the 
DIS spectrum extending from 3650A to 5500A which 
contains Balnier, He I, Mg, and Ca H & K lines (see 
Fig. EJ. The signal-to- noise of this portion of the spec- 
trum ranged from 20 to 60. For both the photometry 
and spectroscopy, we determined the best-fitting mass 
and reddening at a given age using the IDL CURVEFIT 
routine. All fits were done assuming a distance modulus 
of 28.20 (D=4.4 Mpc. lSeth et al.l200fi() . and a metallicity 
of [Fe/H]=-0.4, similar to t he expect ed current gas phase 
metallicity of NGC 4244 ijGarnetti 12002). This is the 
same metallicity derived from composite age fits for a ma- 
jority of thej^ustersjoundby Wal^^^ et al. (2006) . Fur- 
thermorc, Rps sa et all |2006) and Walc her et alJ l|2006f) 
have shown that the inferred age distribution does not 
depend strongly on the metallicity. 

Figure shows the results of our fits for single stel- 
lar populations with ages between 10^ and 10^° years for 
both the spectroscopy and photometry. The top panel 
indicates the of the best-fitting model at each age. 
The spectrum is best fit by populations younger than 
<10^ years, stellar masses of 2-3 x 10^ Mq and decreasing 
extinction with increasing age. The ^ minimum 

at ages of ^0.8-1.0x10^ years (just before the onset of 
the RGB), and then increases rapidly towards older ages. 
This results from the template spectrum becoming red- 
der than the observed spectrum. A shallower minima 
is also seen in the spectra at ~50-100 Myr. 

The fits to the photometric data have two comparable 
minima at ^10* and ~10^ years, with the younger 
age being somewhat preferred. The minima at these two 
ages is more pronounced in the photometric data than 
in the spectroscopic data, perhaps as a result of the pho- 
tometric data's longer wavelength coverage. Combining 
the photometric and spectroscopic results, the minima in 
the joint are at 70 Myr and 0.8 Gyr. The masses from 
the photometric data at these minima are 3-5x10^ M0, 
somewhat more massive than the spectral data. This dif- 
ference may result from light lost outside the 1.5" wide 
slit, uncertainty in the absolute calibration of the spec- 
troscopic data, and/or differences in the assumed IMF. 

Spectral fits from the two joint-x^ minima are shown in 
FigureEl The younger age appears to fit the bluer part of 
the spectrum better, while the older age more accurately 
fits many of the lines and the redder continuum. The 
Ca H line is particularly age sensitive, being much too 
shallow in the younger age and too deep in the older age. 

4.2.2. Multi-Age Fits 

The color difference between the disk and spheroid 
components discussed in M.ll a nd the r ecent results of 
lR.ossa et al.l l|2006>) and .Walcher et all l|2006D strongly 
suggest that NCs contain stellar populations of multi- 
ple ages rather than the single-age populations assumed 
above. We have tried to constrain the ages and masses 
of the multiple stellar populations in NGC 4244's NC 
in two ways: (1) by fitting a weighted average of two 
or three stellar populations of different ages to the spec- 
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Fig. 7. — Photometric and spectroscopic fits of the mass and red- 
dening of the NGC 4244 NC using single-stellar population mod- 
els. Top panel - the reduced values of the best-fitting mass and 
reddening at each age. Middle panel - The best-fitting reddening, 
bottom panel - the best fitting mass. The dashed vertical lines in- 
dicate the two minima in the joint distribution at 70 Myr and 
0.8 Gyr. 



tral and photometric data, or (2) by fitting the spec- 
trum to a co nstant sta r form ation rate (SFR) model 
from Bruzual fc Chariot lj2003|) . These approaches both 
give significantly better fits to the data than single-stellar 
population models, reducing values by 50% or more. 
However, it is also clear that there is no unique solution; 
the data can be fit equally well by many different combina- 
tions of stellar populations and reddening values. Note 
that we assumed a single extinction {Ay) value for all 
components. 

Of the nine NC spectra presented in I Walcher et alJ 
(20Q6i), our spectrum of the NGC 4244 NC most closely 
resembles their spectrum of the NGC 7793 NC. From 
their fits to composite stellar populations, they find that 
the NGC 7793 cluster has an average stellar age of 
~2 Gyr, but has young stars (^100 Myr) that make up 
10% of the mass. This is qualitatively similar to what 
we find below for the NGC 4244 NC. 

Two Age Models: The simplest stellar population 
model consistent with the color maps is a two-age stellar 
population: a young population associated with the disk, 
and an older population associated with the spheroid. 
Taking the two minima seen in the single-stellar popula- 
tion fits at ages of 70 Myr and 0.8 Gyr as a starting point, 
we explored what combination of two populations gives 
the best fit to both the spectroscopic and photometric 
data. For a young component with ages between 50-100 
Myr and an intermediate-age component with ages be- 
tween 0.6 and 1.0 Gyr, the best fits to the spectroscopy 
have Av ^ 0.5, masses of 0.5-3x10^ M0 for the young 
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stellar population, total masses of 2-4x10^ M©, and re- 
duced values of ~6. However, for the photometry, 
the best fits favor an Ay > 1 and masses dominated by 
the young component. These photometric fits clearly do 
not reproduce the shape or lines in the spectra, thus we 
restricted ourselves to lower Ay solutions. When Ay is 
restricted to be less than one, the photometric fits are 
in relatively good agreement with the spectroscopic fits. 
However, the implied stellar masses for both the young 
and intermediate-age components are up to two times 
larger than the spectroscopic best-fit models, similar to 
what was seen in the single-age fits. Overall, we find that 
ages of 0.1 and 1.0 Gyr do the best job of matching both 
sets of data. The spectroscopic best fit at these ages is 
shown in Figure El and has Ay = 0.46, 1.9x10^ M© of 
0.1 Gyr old stars and 3.1x10^ Mq of 1 Gyr old stars. 

Based on the color map of the nuclear cluster, we ex- 
pect the younger stellar component in the two-age fits 
be associated with the disk. We can test the consis- 
tency of this idea by comparing our spectral fits to the 
luminosities of the disk component decomposed from the 
HST data in Assuming that the disk is dominated 
by '-TO Myr old stars, the expect ed F606W-F8 14W color 
from the Padova models is 0.25 l)Girardill2006j) . The ob- 
served disk component color is 0.4, requiring an extinc- 
tion Ay = 0.47 mags. Using this extinction and mass- 
to- light ratios from the Padova models, we derive a mass 
for the disk component of 1.2x10^ M©. Both the mass 
and extinction derived from the luminosities of the disk 
component in the HST data agree with the young stellar 
mass and extinction determined in the two-age spectral 
fits. This agreement strongly suggests a link between 
the blue disk component and the young stellar compo- 
nent that appears to be required to fit the blue spectral 
features. 

Three Age Models: We also explored the sensitivity 
of our fits to the inclusion of a third, old (10 Gyr) stellar 
population. We find that a large mass of old stars is 
consistent with both the photometric and spectroscopic 
data. Best-fits to both data sets have low Ay (<0.4), 
total masses of '-^lO^ Mq with 90% of that mass in the 
old component, and '-~^10^ M© in the young component, 
similar to the two-age fits. The inclusion of the older 
component significantly improves the photometric fit to 
the data (reduced = 1-7), suggesting older stars may 
indeed be present. However, the spectral fit is insensitive 
to this old component; the best-fit spectrum is nearly 
identical to the two-age fits, despite 90% of the mass 
being contained in the old component. 

The effects of disk contamination on the photome- 
try and spectra should be minimal given the promi- 
nence of the NGC 4244 NC. Futhermor e, by compar- 
ing HS T/STIS and ground-based spectra. FWalcher et al.l 
l)200fi|) show that the stellar population immediately sur- 
rounding the nuclear cluster in late-type galaxies seems 
to have a p o pulati on similar to the cluster itself. Both 
lR,ossa et all ^06) and 'W alcher etal] (|2006) find some 
evidence for an underlying old stellar component in many 
of their nuclear cluster spectra. We note that all the 
masses quoted here refer to the initial mass in stars re- 
quired to fit the current spectral data. Due to stellar 
evolution, the current mass of a cluster dominated by 
old stars would be '^2x lower than the initial mass, given 
reasonable initial mass functions. 



Constant SFR model: Using the iBruzual fc Chariot 
l|200,'D codes, we generated a model spectrum assuming a 
constant SFR from 12 Gyr to the present. The best-fit of 
this model to the data gives a current mass of the cluster 
of 6.2x10^ Mq, and a reddening of Ay = 0.72. The re- 
duced of the fit is 5.7, within 10% of the reduced of 
the best fits obtained for the two- and three-age models. 
The spectral shape of the constant SFR model appears to 
be somewhat redder than the data, suggesting either the 
need for relatively more young stars than in the constant 
SFR model, or for a more realistic metallicity evolution, 
in which the older stars are bluer and more metal-poor . 
The former point is supported bv I Walcher et al.l l)2006D . 
who have shown that in their sample of clusters, the star 
formation rate in the past 100 Myr is an order of mag- 
nitude greater than the star formation rate between 0.3 
and 20 Gyr. However, this enhancement of recent star 
formation may result from a bias in their sample towards 
brighter nuclear clusters. 

4.2.3. Summary of Stellar Populations in the 
NGC 4244 NC 

We draw the following conclusions from the spectral 
and photometric fits of the NGC 4244 NC: 

1. Multiple stellar populations definitely provide im- 
proved fits over single-stellar populations for both 
the spectroscopy and photometry of the NGC 4244 
NC. However, many different combinations of 
masses and ages fit the data well. 

2. The best-fitting single stellar populations have ages 
of either --70 Myr or -0.8 Gyr. 

3. The total current mass of the cluster appears to be 
between 0.2 and IxlO^M g. 

4. If we assume only two ages are present, the 
youngest (< 100 Myr) stellar component appears 
to have ~10^ Mq of stars. This mass is consistent 
with the luminosity of the disk component decom- 
posed using the morphological fits in ^ 

5. For all multi-age fits, the mass of stars younger 
than 100 Myr is a small fraction (< 5%) of the 
total mass of the cluster. 

6. Stars older than 1 Gyr may be present in large 
numbers (—10^ Mq). Inclusion of an old stellar 
population particularly improves our fits to the 
SDSS and 2MASS photometric data. 

7. A constant SFR fits the data as well as assuming 
two or three discrete ages. 

5. A DYNAMICAL MASS ESTIMATE FOR THE NGC 4244 
CLUSTER 

Emission lines, including Ha, [Nil], [SII], and [Oil] 
are seen in our spectrum of the NGC 4244 NC. For the 
Ha line, the emission component is redshifted relative to 
the underlying stellar absorption lines (see Fig. ^ . We 
utilize this offset between the emission and absorption 
components to obtain a dynamical mass estimate for the 
NGC 4244 cluster. 

Examination of the color map for NGC 4244 in Fig- 
ure [S] shows a very blue compact source located along 
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Fig. 8. — Spectral lines used to derive the relative offset between 
emission and absorption. Dark lines show the actual spectrum, 
while lighter lines show the best fitting gaussian or double-gaussian 
model. The number in the title of each panel gives the fitted ve- 
locity and error bar. The spectra has a S/N~35/pixel. 



the NC plane on the left (NE) side, 0.88"(18.6 pc pro- 
jected) from the NC center. This compact blue object 
has an F606W-F814W color of -0.3, which is -0.4-0.6 
mags bluer than the disk and spheroid components of 
the cluster. Assuming moderate reddening and blending 
with the outer isophotes of the cluster, the color of the 
blue knot is consistent with this object being an HII re- 
gion. HII regions in our filter combination appear very 
blue because there are strong emission lines located in 
the F606W filter and not in the F814W filter. Measure- 
ments of HII regions in NGC 55 (the nearest galaxy in 
our sample) suggests HII regions have typical colors of 
F606W-F814W=-0.2 with a -0.2 mag scatter. 

This likely HII region falls within the slit used for our 
nuclear cluster spectrum. Therefore, we identified this 
object as the likely source of the emission lines in our 
spectra. We verified this by extracting spectra from the 
left (NE) and right (SW) sides of the cluster. The left 
side spectrum shows significantly enhanced emission lines 
relative to the right side. However, because the seeing 
was —1", we were unable to spatially separate the emis- 
sion component from the nuclear cluster. We now pro- 
ceed under the assumption that the emission lines in our 
spectrum are associated with this compact blue object. 

To determine the relative velocity of the emission and 
absorption components we fit the [SII] doublet and [Nil] 
line to determine the emission line velocity, and fit the 
Nal doublet to determine the absorption line velocity. 
We obtained a weighted mean velocity of 270.5±5.4 
km/sec for the emission lines and 246.3±4.3 km/sec for 
the Nal absorption line assuming air wavelengths for 
the lines from the Atomic Line List^. This absorp- 
tion velocity match es the system ic velocity of NGC 4244 
(244.4±0.3 km/sec, '011inE"1996") to within the error bars, 
as would be expected for the nuclear cluster. We note 
that this HII region is counter-rotating with respect to 

^ http:/ /www. pa. uky.cdu/~ peter/atomic/ 



the galaxy disk. 

FigurelHlshows the fitted emission and absorption lines, 
as well as the best fit of absorption -I- emission for the 
Ha line after fixing the velocities of the absorption and 
emission to the values derived from the Nal, [SII], and 
[Nil] lines. This Ha line fit does a good job of replicat- 
ing the observed line shape thereby suggesting that the 
velocities we have derived are correct. 

Assuming the projected radius of the HII region is close 
to the true radius, and that it is bound to the NGC 4244 
NC nuclear cluster, the Av — 24.1 ± 7.0 km/sec offset 
between emission and absorption implies a nuclear clus- 
ter mass of 2.5l]^'2 ^ lO^M© within a radius of 19 pc. If 
the object is indeed rotating around the cluster, but is in 
front of or behind the cluster (thus making its projected 
distance larger), then the mass we derived is a lower limit 
on the mass of the cluster. 

This rough estimate of the dynamical mass of 
NGC 4244's NC is consistent with the stellar masses de- 
termined for the one- and two-age fits in M.2\ It is also 
consistent with the constant SFR or three-age model fits 
(M— 6x10^ Mq) if the HII region is seen in projection 
and lies —25 pc from the center of the nuclear cluster. 

The estimat ed mass is also simila r to the typical mass 
determined by'Walcher et al.' (2005) for late- type galaxy 
nuclei. They derived masses for nine nuclear clusters 
from the iBoker et alJ p002il sample by combining size 
and velocity dispersion measurements. Their total range 
of masses is 8x10^ M0 to 6x10'' Mq, with seven of the 
nine clusters having masses between 0.1 and 1x10^ M0. 
We note that they have determined these masses as- 
suming spherical symmetry; if these clusters have disk 
components similar to the ones we have ob served here, 
the masses derived bv lWalcher et all lj2005() may be un- 
derestimates, since they do not take into account the 
rotational support that would be expected in a disk. 
Given that their clusters c ontain young stellar popula- 
tions ijWalcher et al.l 12006(1 . this effect could be signifi- 
cant. 

6. NGC 4206: POSSIBLE INDICATION OF AN AGN 
COMPONENT 

In this section, we present evidence that the nuclear 
cluster in NGC 4206 may contain an AGN component. 
This evidence is very tentative and requires confirmation 
with follow-up observations. 

An SDSS spectrum (MJD 53149, Plate 1612, Fiber 
603) exists for the central regions of NGC 4206. Due 
to the 3" diameter (—250 pc) aperture, this spectrum is 
dominated by the light from the central regions of the 
galaxy and not by the nuclear cluster. In addition to 
a dominant old stellar spectrum, narrow Ha, [Nil], and 
[SII] emission are clearly seen. No broad emission lines 
are seen. Examining the HST color map, the center of 
the nuclear cluster is the bluest area anywhere within the 
SDSS spectral aperture. Since emission line sources are 
very blue in our filter combination, the central portion 
of the nuclear cluster is the likely source of the emission 
lines. The emission lines are not offset in velocity from 
the absorption spectrum (as was seen in NGC 4244) fur- 
ther supporting their nuclear origin. 

The ionizing source of the observed emission lines could 
be either star formation or AGN activity at the center 
of the nuclear cluster. These two ionizing sources can be 
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distinguished using the hne ratio of [Nil] and Ha. Be- 
cause the Ha line is contaminated by stellar absorption, 
we model the underlying stellar continuum with a single 
stellar population '^3 Gyr in age with Ay ~ 0.5. After 
subtracting this spectrum, the line ratio of the remain- 
ing emission lines is log([NII]/Ha)~ —0.3. Based on th e 
recent analysis of SDSS spectra by ijObric et al.ll2006() . 
this line ratio suggests that the source is a "mixed" star- 
burst/ AGN. NGC 4206 may therefore host a small AGN 
within the core of the nuclear star cluster. If confirmed 
by follow-up observations, this would make it only the 
second late-type spiral known to have both a SMBH and 
massive nuclear star cluster, the other being NGC 4395, 
an Sd galaxy with Mn ^ —17.5 hosting a n uclear cluster 
with Mi -10.1 ijMatthews et al.l llQQQl. and a black 
hole with a mass of -lO^ - 10^ Mq (Filippcnk o fc Hoi 
|2003|) . In earlier type spirals, other examples exist of 
the co-existen ce of nuclear star clusters and AGN (e.g. 
iThatte et all irQQ?: Davics ct al. 2006). We plan to test 
for the presence of AGN components in all of our multi- 
component nuclear cluster candidates with near-IR inte- 
gral field unit spectroscopy using adaptive optics. 

7. DISCUSSION 

We now discuss how our data constrains the formation 
mechanisms of nuclear clusters. We suggest a forma- 
tion mechanism in which stars in nuclear clusters form 
episodically in disks and then over time lose angular mo- 
mentum or are heated vertically, ending up in a more 
spheroidal distribution. After discussing this mechanism 
in some detail, we examine the relation of nuclear clus- 
ters to other massive star clusters and supermassive black 
holes. 

7.1. Formation Mechanisms 

Our observations of multi-component clusters place 
significant constraints on nuclear cluster formation mech- 
anisms. The presence of young NC disks (<1 Gyr) 
aligned with the disks of the host galaxies strongly sug- 
gests that they are formed in situ from gas accreted 
into the nuclear regions. We have shown in SjHl2 that 
these disks may be a common feature of nuclear clus- 
ters. Such disks cannot be formed via the accretion of 
globular clusters by dynamical friction. While globu- 
lar cluster accretion very likely occurs, this mechanism 
has bee n shown to b e relatively inefficient in late-type 
galaxies l|Milosavlie vic 2004). We propose here a scenario 
in which the dominant formation mechanism for nuclear 
clusters is via in situ formation in disks. Although the 
mechanism will need to be tested with detailed simula- 
tions, it is consistent not just with our data, but with 
both the evidence for young and co mposite stellar popu - 
lations in many nuclear clusters fe.g lWalcher et all2006j) . 
and the direct detection of a molecula r gas disk coinci- 
dent with the nuclear cluster in IC 342 ijSchinnerer et al.l 
11)03). 

On a large scale (~1 kpc), both dissipative merging 
of molecular clouds and magnetorotational instabilities 
have been shown to be capable of supplying the nu- 
cleus of a late-tTOC galaxy with sufficient gas to form 
nuclear clusters (|Milosavlievic 2004; Bckki et al. 20061). 
Our proposed in situ formation mechanism is concerned 
with what occurs once the gas reaches the center. Our 
model is also fully consistent with the suggestion by 



iWalcher et alJ 1)200611 that an initial seed star cluster is 
required for subsequent nuclear star cluster formation. 

The details of an in situ formation mechanism can be 
constrained by our observations. To explain presence of 
redder/older spheroids in the multi-component clusters 
and the absence of disk components in a majority of our 
nuclear clusters, we propose a model in which nuclear 
cluster disks form episodically and are transformed into a 
rounder distribution over time due to loss of angular mo- 
mentum or dynamical heating. There are two important 
timescales in this model: (1) the period of time between 
star formation episodes, and (2) the time it takes stars in 
a disk morphology to end up in a spheroidal distribution. 
We examine these two timescales in detail below. 

7.1.1. Period Between Star Formation Episodes 

A number of authors have previously suggested that 
star formation in the nuclear regions of late-type spirals 
may be episodic. The physical motivation for episodic 
star formation is that a massive star formation event 
will feedback on the local environment preventing star 
formation for the next ~10® years (|Milosavlicvic 2004). 
Although this feedback may destroy the gas disk and 
shut off star formation, it does not disrupt the disk of 
recently formed stars, whose lifetime we discuss below in 
U\1.2. 

Assuming that the present epoch is not a special time 
in the life of the nuclear clusters, one can estimate the 
peri od on which star for r nation recurs. Observation- 
ally, 'Davi dge fc Courteaul 1)2002)) argue for a period of 
^0.3 Gyr based on t he detection frequency o f Paa emis- 
sion in galaxy nuclei. iSchinnerer et alJ l)2003)l also deter- 
mine a recurrence period of a few hundred Myr to 1 Gyr 
based on the gas accretion rates observed in IC 342 's nu- 
cleus. More recently, from analysis of the youngest stellar 
components in the spectrum, Walcher ct al. ( 2006) has 
found the typical time between star formation bursts to 
be ~70 Myr with a typical mass of roughly 1.6 x 10^ Mq. 

Our observations of the NGC 4244 nuclear cluster are 
consistent with these previous estimates. Based on the 
^10^ Mq of stars in the young component f ti4.2|l . ap- 
proximately 25 episodes of disk formation must have oc- 
curred to build up the 2.5x10^ Mq of the cluster (©• If 
this occurs over a Hubble time, then the period between 
bursts is ^0.5 Gyr. 

7.1.2. Timescale of Disk Transformation 

Given that nuclear cluster disks are observed in <50% 
of nuclear clusters, and that these disks all have stel- 
lar populations younger than a Gyr, a mechanism must 
exist that destroys or disrupts this disk on a timescale 
shorter than ^1 Gyr. We suggest that angular mo- 
mentum loss or dynamical heating causes the disks 
to be transformed over time, thus creating the more 
spheroidal, older components seen in our nuclear clus- 
ters. Although mechanisms of angular momentum loss 
have not been studied previously in nuclear clusters, 
these objects are dynamically similar to globular clus- 
ters, for which the theory of rot ation and rotational 
evolution has been wcU-dcvcloped fAgekianlll958t iKind 
0^61; Shaoiro & Marchant 1976: Fall fc Fren^ ll9pF 
iDavoust fc Prugniel 1990; Longa retti fc Lagoutel ll9971 

As they age, star clusters flattened due to rotation 
will become increasingly round due to the preferential 
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evapo ration of high angular momentum stars l|AgekiarJ 
Il958j) . We can compare the timescale of this angu- 
lar momentum loss to the expected lifetimes of our 
disks. The timescale over which the ellipticity of a rotat- 
ing stellar clust er decreases is > 40x the median relax- 
ation time, Trh l|ShaDiro fc Marehant 1976; Fall fc Re^ 
Il985() . For globular clusters in the Milky Way, Trh is 
typically ~10* years. However, this timescale is longer 
for more massive clusters lik e lo Cen, with r^h = 3.5 Gyr 
ijShapiro fc Ma rehant 1976^. For the NGC 4244 NC we 
have calculated a median rel axation time for the whol e 
cluster of -2 Gyr (Eq. 8-72 lBinnev fc Tremain3ll987t) . 
For just the disk component, assuming it has a mass 
of ^10^ M0 and a half mass radius of ~3 pc, we get 
Trh = 4 X 10^ years. These relaxation times suggest that 
the timescale for evaporation of angular momentum in 
NGC 4244 is in the range of 16-80 Gyr. This timescale 
is significantly longer than the 100s of Myrs over which 
we expect the disk to be transformed into a spheroidal 
component. 

It therefore appears that angular momentum loss due 
to processes internal to the clusters are unlikely to trans- 
form the disk components into spheroids on sufficiently 
short timescales. Instead, we must look to external pro- 
cesses that can affect the the dynamical evolution of the 
nuclear cluster. We have identified a few alternative 
mechanisms for transforming the stellar disk component 
into a spheroid. These include: (1) small angular off- 
sets of the nuclear disk from the plane of the galaxy, and 
(2) slight positional offsets of the cluster from the center 
of the galaxy potential, (3) torques on the stellar disk 
from successive gas accretion. A forthcoming paper will 
discuss this problem in more detail. 

7.1.3. Rings 

We have been assuming thus far that stars are formed 
initially in disks, as suggested by our observations of 
NGC 4244 and NGC 4206. However, in 101.3, we showed 
that one of the three multi-component NCs, in IC 5052, 
was best fit by a spheroid -I- ring model with inner and 
outer ring radii of 6.8 and 11.8 pc. We discuss here how 
such a ring could form and examine rings in context of 
our proposed nuclear cluster formation model. 

Base d on an elongated near-infrared component in 
IC 342, iBoker et all i)1997|) suggested that the observed 
molecular ring (with radius^ of ^65 pc) results from an 
inner Lindblad resonances (ILR) with a stellar bar. If 
a similar mechanism is at work in forming the ring in 
IC 5052, it would require the star cluster itself to have a 
barred disk morphology. From HST imaging of nuclear 
clusters in fa ce-on spirals, there is no evidence for bars on 
these scales l|Boker et al.ll2004j) . A ring could also result 
from tidal stripping of an infalling molecular cloud or 
star cluster. Assuming a pre-existing cluster with mass 
of 3 X 10^ Mq and a 10^ M0 molecular cloud, the molec- 
ular cloud would be tidally stripped between radii of 8 
and 20 pc for molecular cloud number densities ranging 
between 5 x 10"* and 5 x 10'^ cm""^. Thus, tidal stripping 
of a molecular cloud is a feasible mechanism for creating 
a ring such as the one in IC 5052. 

Episodic star formation occuring at least occasionally 

^ The radius of their ring was recalculat ed assuming the Cepheid 
distance to IC 342 of 3.3 Mpc derived bv lSaha et'aTT<200^ . 



in rings is consistent with the proposed formation model. 
The ring in IC 5052 is internal to the cluster itself, so 
after loss of angular momentum or dynamical heating, 
its stars would still contribute to the stellar distribution 
of the spheroid. The radius of the ring depends on the 
current mass of the star cluster in both mechanisms for 
ring formation described above. This correlation might 
provide a natural explanation for the clear correlation 
of nuclear cluster si ze and luminosity observed in dwarf 
elliptical galaxies bv lCote et al.l |2006). A similar expla- 
nation could work for formation in disks if the size of the 
disks depends on the mass of the nuclear cluster. De- 
tailed modelling of disk/ring formation and dissipation 
would be required to test the feasibility of this idea. 

7.1.4. Other Formation Scenarios 

Our proposed scenario naturally explains the forma- 
tion of multi-component clusters. However other scenar- 
ios may also be consistent with our observations. While 
our observations seem to require some in situ formation 
to explain the presence of young stellar disks, they do 
not rule out the possibility that the multi-component 
nuclear clusters were f ormed in part from the merging of 
globular clusters (e.g. iTremaine et alJll97,'Tl : iLotz et all 
I2OOIID . In this "hybrid" scenario, the young disk com- 
ponents may be a late-time veneer on an already mas- 
sive cluster formed from merged g l obula r clusters. How- 
ever, we note that iMilosavlievicI l)2004D has calculated 
that the timescale for dynamical friction in late-type spi- 
rals is > 3 X 10^ years even for relatively massive clusters 
(10^ M0) located within the central kpc of the galaxy. 
This quantity increases for less massive and more distant 
clusters. He concludes that this timescale is thus too long 
to be the dominant mechanism for nuclear cluster forma- 
tion, and shows that magnetorotational instabilities can 
transport sufficient quantities of gas into the galaxy cen- 
ter to account for the observed NCs in late-type spirals. 

7.2. Relationship to Other Massive Star Clusters 

Our observations have implications beyond just the 
formation of nuclear star clusters in late-type galaxies. 
We briefiy review here the possible connections between 
late-type spiral nuclei and dwarf elliptical nuclei, and the 
evidence that some massive globular clusters and ultra- 
compact dwarfs are stripped nuclear clusters. 

7.2.1. The relation between spiral and dE nuclei 

The recent study bv lCote et al.l 1)20061) shows that the 
nuclei of dE and late-type spiral galaxies have compara- 
ble sizes, luminosities, and frequencies of detection (see 
also Fig. 121). They argue that this similarity suggests a 
nuclear cluster formation mechanism that does not de- 
pend strongly on galaxy properties. However, the forma- 
tion scenario suggested by our observations for nuclear 
clusters in late-type galaxies cannot take place in the gas- 
free dE galaxies. We therefore interpret the similarity 
of the nuclei in these galaxies as support for dE galaxies 
evolving from low-mass s piral and irregular galaxies (e.g. 
iDavies fc PhilliDDsll988D . In this scenario, a dE,N galaxy 
would result from a nucleated late-type spiral galaxy af- 
ter depletion of gas through s tar formation or galaxy har- 
rassment ( Moore et al. 1996). If so, then the ages of the 
nuclear clusters in dE galaxies can constrain the epoch 
when these systems were stripped. 
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7.2.2. Massive star clusters as stripped galaxy nuclei 

There is direct evidence that some massive globular 
clusters are associated with stripped galaxies. The galac- 
tic globular clu ster M54 appears to be th e nucleus of 
the Sag dSph l)Lavden fc Saraiedinil l2000j) . and six of 
the most massive globular clusters in NGC 5128 appear 
to be surrounded by extratidal light (Harris et al. 200^ 
IMartini fc Ho"2004'l. 

Stripped nuclear clusters have also been invoked to 
explain the multiple stellar populations o bserved in the 
most massive nearby globular clusters fe.g lGnedin et all 
l2002j) . These observations include the large metal- 
licity spre ad ( Sniith| I2004f) and double main-sequence 
ijBedin et al 2004) in co Cen, and the spr ead in red gi- 
ant b ranch colors for Gl in Andromeda ijMevlan et al.l 
1200 ID . More recently, iHarris et all l)2006j) have observed 
a trend towards redder colors with increasing luminosity 
for the brightest blue globular clusters in massive ellip- 
tical galaxies. This trend is interpreted as evidence for 
self-enrichment in these clusters pointing to an extended 
star formation history more typical of nuclear clusters 
than globular clu sters. A similar r edward trend is seen 
for dE nuclei bv ICote et alJ l)2006f) . Self-enrichment is 
naturally explained in an in situ formation mechanism, 
but is not easily explained if these clusters are formed by 
merging less massive clusters via dynamical friction. 

The ultr acompact dwarf galaxies recently found in 
the Fornax ( PhillioDS ct al. 2001: Dri nkwater e t al.'2003') 
and Virgo ( Jon es_et_aL 2006) clusters may also be re- 
lated to nuclear star clusters. One possible mechanism 
for their formation is "galaxy threshing", in which nu- 
cleated galaxies lose t heir envelopes leav ing behind only 
their nuclei as UCDs ( Bckki et al. 2001). In the Mi vs. 
log(re//) diagram in FigurcEl NGC 4206's NC has suffi- 
cient size and luminosity, even after modest age fading, to 
be a UCD progenitor. The galaxy threshing sc enario has 
also been shown to be feasible bvlCote et al.l (2006) for 
dwarf elliptical nuclei, although iDe Prooris et al.l (l|200fi^ 
suggests that dE,N are in general more compact than 
UCDs. One of the galaxy nuclei they study does have 
similar properties to the UCDs, as well as possessing spi- 
ral arms which indicate the galaxy may be an anemic 
spiral. 

7.3. Connection to Black Holes 

As residents of the same neighborhoods, it is natural 
to suspect that there may be some connection between 
nuclear clusters and supermassive black holes (SMBHs). 
The presence of SMBHs in the nuclei of >30 galax- 
ies have been confirmed through dy namical measure- 
ment s of their masses (reviewed by iFerrarese fc ForJ 
l2?ffl5l) . The masses of SMBHs are known to correlate 
with the mass and/or velocity dis persion of the bulge 
component (t he M, — a re lation, IFerrarese fc MerrittI 
iOOO; GcbhardTeEaD MM)- Recent work on nu- 
clear clusters in both spiral and ellip tical galaxie s have 
shown a remarkably simila r trend i Rossa et al.l 120061: 
Coteet a] .II20M IFerrarese et al. 2006HWehuer fc Ha,rnFi 



2006|). For a large sample of nucleated elliptical galax- 
ies. ICot^ t al. (2003) have shown that the nucleus ac- 
counts for 0.30±0.04% of the total galaxy mass. This 
fraction is very simila r to the fraction of g alaxy mass 
found in the SMBHs (jFerrarese et al.ll2006(l . For this 



reason. iCote et al.[ (I2006D . IFerrarese et al.l (|2006D and 
iWehner fc HarrisI l)2006l) siiggest that nuclear clusters 
and SMBHs be thought of as a single-class of Central 
Massive Objects. They suggest a transition at ^^10^ 
Mq between nuclear clusters and SMBHs. However, the 
detection of numerous lower mass SMBHs at the cen- 
ters of early and late-ty pe dwarf galaxies (jBarth et alJ 
I2OO5H G reene et al.l2005ll . and the co-existence of nuclear 
star c lusters and SMBHs in NGC 4395 (Filippcnko fc Hd 
12003(1 and perhaps in NGC 4206 ([0, suggest that this 
transition may not be abrupt. 

The formation of black holes in the centers of nu- 
clear clusters may be related to the intermediate mass 
black holes (IMBHs) found in the centers of massive non- 
nuclear star cl usters. A 1.7x lO^M© black hole has been 
reported by .Geb hardt et al.l ((20051 ) in the globular clus- 
ter Gl in Andromeda (which is also a possible stripped 
nuclei), based on stellar dynamics. A lO^-lO'' Mq black 
hole has also been identified in one of M8 2's young mas- 
sive c lusters, based on X-ray observations ('Pa truno et alJ 
12006(1 . These black holes are thought to form either by 
early merging of massive young stars within the clus- 
ter, or by dynamical interaction in the core of the clus- 
ters over long time periods (Mil ler fc Colber3l2004l and 
references therein). If IMBH formation is common in 
dense massive clusters, then nuclear clusters should fre- 
quently host black holes as well. Subsequent mergers of 
many IMBHs could then for m sup ermassive black holes 
(SMBHs) ((Miller fc ColbertI I200I . As the mass of the 
SMBH increases, these mergers will likely destroy the nu- 
clear clusters (see §5.2.2 of Cote et al..2006,) . distributing 
their stars into the bulge. 

8. SUMMARY 

We have examined the nuclear regions of fourteen 
nearby galaxies using HST/ACS images in the F606W 
and F814W filters. We find nuclear cluster candidates 
within 2" of the photocenter of nine of these galaxies. We 
have also obtained multi-band photometry and a spec- 
trum of the NGC 4244 nuclear cluster, the nearest object 
in our sample. Analyzing these nuclear clusters, we find 
the following: 

1. From the edge-on perspective, three of the 
nine cluster candidates (IC 5052, NGC 4206, 
NGC 4244) have morphologies suggesting they 
have both disk and spheroid components. These 
components are very compact with scale lengths 
and half-light radii of 2-30 pc. Star clusters with 
multiple visible components have not previously 
been seen. 

2. The three multi-component nuclear clusters have 
elongations oriented to within ~10° of the major- 
axis of the galaxy. 

3. The magnitudes and sizes of the clusters in our 
edge-on sample are very similar to th ose observed 
in ot her face-on, late-type galaxies ("B oker et alJ 
|2004() and dE galaxies (Cote et al. 2006). 

4. The spheroids of the multi-component nuclear clus- 
ters are 0.3-0.6 magnitudes redder than the disks 
in F606W-F814W. This color difference implies the 
disks are populated with stars <1 Gyr in age. 



Nuclear Star Clusters in Edge-on Spirals 



17 



5. Fitting single-stellar population models to the 
spectral and photometric data of the NGC 4244 
nuclear cluster, we find it is best-fit by ages of 
'~70 Myr or ~0.8 Gyr. The implied total mass 
of the cluster for is ~3xl0^ M©. 

6. A combination of two or more stellar populations 
significantly improves the fit to the NGC 4244 nu- 
clear cluster spectral and photometric data. In the 
multi-age fits, the youngest (< 100 Myr) stellar 
population matches the luminosity of the disk com- 
ponent, and appears to make up ^5% of the total 
mass of the cluster. As much as ^10^ Mq of old 
(10 Gyr) stars may also be present in the cluster. 
A constant star formation rate over the last 12 Gyr 
also provides a surprisingly good fit. 

7. Using an HII region 19 pc from the center of the 
NGC 4244 NC, we measure the velocity offset of 
emission and absorption components in our spec- 
trum to obtain a lower limit to the dynamical mass 

-t-1 7 

for the nuclear cluster of 2.5_j2 x lO^M©. 

8. Analysis of emission lines in the NGC 4206 nucleus 
suggest that it may host an AGN as well as a stellar 



cluster. 



These observational results strongly support an in situ 
formation model for nuclear star clusters. We suggest 
a specific model in which stars are formed episodically 
when gas accretes into a compact nuclear disk. After 
formation the stars gradually lose angular momentum 
and eventually end up in a spheroidal component. We 
show that a duty cycle for star formation of ~0.5 Gyr 
is consistent with our and previous observations. The 
angular momentum must also be lost on this timescale, 
which we suggest may occur due to interactions of the 
nuclear cluster with the surrounding galaxy potential. 
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